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The effect of polydispersity on the freezing transition of 
hard spheres is examined within a moment description. At 
low polydispersities a single fluid-to-crystal transition is re- 
covered. With increasing polydispersity we find a density 
above which the crystal melts back into an amorphous phase. 
The range of densities over which the crystalline phase is sta- 
ble shrinks with increasing polydispersity until, at a certain 
level of polydispersity, the crystal disappears completely from 
the equilibrium phase diagram. The two transitions converge 
to a single point which we identify as the polydisperse ana- 
logue of a point of equal concentration. At this point, the 
freezing transition is continuous in a thermodynamic sense. 

PACS numbers: 05.70.Fh, 64.60.Cn, 64.70.Dv, 82.70.-y 



Freezing and melting are probably the most common 
and striking physical changes observed in everyday life. 
All experiments, to date, demonstrate that the crystalli- 
sation of a simple liquid is a first-order transition, in three 
dimensions. So for instance, the sharp Bragg peaks of the 
crystal, which reflect the long-range spatial modulation 
of the density p(r) and which distinguish a crystal from 
a liquid, disappear abruptly as a crystal melts Q. This 
sharp microstructural change is also mirrored by discon- 
tinuities in the first derivative of the free energy so that 
experimentally, melting is accompanied by a finite den- 
sity and entropy change. 

Although the experimental situation is clear, in an 
early analysis Landau argued that, under certain 
conditions, a crystal can transform continuously into a 
liquid. In a simple Landau-Alexander-McTague theory 
|pf the excess free energy of the crystal (relative to the 
isotropic liquid) has the following form: 
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where the order parameters ng are the Fourier compo- 
nents of the crystal density, p s {f) — Ps + Sp{r), at the 
reciprocal lattice vector G (p s is the uniform crystal den- 
sity) and the coefficients of the expansion are analytic 
functions of the temperature T and pressure P. Eq. |l| 
contains cubic terms because the order parameter sets 
{no} and {—no} describe physically distinct crystals 



with different energies. As a consequence the freezing 
transition is generally first-order. However since both T 
and P can be independently varied the possibility exists 
that r and 1*3 can be made to vanish at a single point in 
the T—P plane. At the resulting Landau point the liquid- 
solid transition is continuous in a mean-field description 
Q4J . Landau theory makes two further distinctive predic- 
tions. First the Landau point must lie at the intersection 
of, at least, three first-order lines of transitions || which 
separate the liquid from two conjugate crystalline phases, 
C+ and C_, with identical symmetry but which differ in 
the sign of Sp(r). Second, in three dimensions, symmetry 
considerations should uniquely favour a bcc structure [p| . 
In spite of these interesting predictions it is not clear if, 
in a liquid-solid system, the point at which the cubic co- 
efficient U3 vanishes is experimentally accessible. On the 
face of it, one of the most promising candidates is a sys- 
tem of polydisperse hard spheres where the constituent 
particles have different sizes. The freezing of polydis- 
perse hard spheres has been studied extensively in recent 
years pH14[ motivated, in part, because it is a realis- 
tic model of a colloidal suspension Ea|. These studies 
have focused mainly on the effect of size polydispersity 
a, defined as the ratio of the standard deviation to the 
mean of the diameter distribution, upon the fluid-solid 
transition. Calculations have been made using a variety 
of theoretical and computational techniques, for various 
size distributions, and in both two and three dimensions. 
Yet the picture that has emerged is remarkably similar. 
On increasing a, from zero the density discontinuity at 
the transition Ap = p s — pi decreases, vanishing alto- 
gether at a "terminal" polydispersity Q , a = a t , above 
which no liquid-solid transition is found. A number of 
key questions have however been left unanswered. First, 
why do the densities of the coexisting phases converge as 
a — ► <7t? If the liquid-solid transition is continuous then 
the singularity at at must correspond to a Landau point. 
The phase diagram should therefore contain two crys- 
tal phases, in contradiction with the theoretical work to 
date. Furthermore while the C+ crystal has the normal 
bcc structure with spheres at the cube corners and cen- 
tre, the C_ crystal has particles at interstitial sites. The 
unfavourably low packing of the C_ crystal ((p m ~ 0.20) 
makes it unlikely that this phase could be important in 
a dense system. If the vanishing density discontinuity at 
at is not critical in origin, then what is its true nature? 



And finally, why is the polydisperse phase behavior ap- 
parently universal? In this letter we reexamine the freez- 
ing of polydisperse hard spheres using simple mean-field 
models for the polydisperse crystal and liquid phases. 
Our results suggest that the polydisperse solid-liquid tran- 
sition at at is not critical. We show that the vanishing of 
the density discontinuity at the terminal polydispersity 
is a consequence of a re-entrant solid-liquid transition in 
a polydisperse system. 

Our model consists of N hard sphere particles in a 
volume V, at an overall density of p = N/V. Each 
particle has a diameter R drawn from a distribution 
p{R) so that p — f dRp(R). The distribution p(R) is 
conveniently characterised by the set of generalised mo- 
ments mi = J dRp(R)u>i(R) where the weight function 
Wi(R) = (R/R—iy. The zeroth moment is simply the to- 
tal number density p. The "shape" of the diameter distri- 
bution, p(R) = p(R)/p, is taken here for simplicity, as the 
Schultz distribution, p(R) =_ ^ a R a - 1 ocp{-~fR)/Y{a) 
with a = 1/c 2 and 7 = a/R. The (excess) chemical 
potential p ex (R) in a polydisperse system is in general a 
complex and unknown function of the particle size. But 
with the assumption that there is no critical point at 
at the excess chemical potential must, first of all, be an 
analytic function of R. Formally, p ex (R) may be calcu- 
lated from the probability, W(R) : for insertion |17J of a 
test sphere of diameter R. At large R, the leading term 
in p ex (R) is the PV work required to generate a cavity 
sufficiently large to accommodate the test sphere. This 
contribution varies as R 3 . Motivated by this we assume 
that in a hard-sphere crystal or fluid p ex (R) has the sim- 
ple analytic form 



p ex {R) = -k B T\nW{R) 
« Aq + \\R + A2-R 



A3.R 3 . 



(2) 



where consistency demands that the coefficients Aj de- 
pend only on the four moments too, . . . TO3 of the poly- 
disperse distribution [|12| . Two of the four unknown co- 
efficients may be determined from the known small and 
large R limits of W(R). This fixes /3A = - ln(l - <f>) and 
A3 = ?P with <fi the volume fraction and j3 = 1/ksT. 

Having specified the general form expected for p ex (R) , 
we now outline the calculation of the size-dependent 
chemical potential in the crystal. From Eq. || the prob- 
ability to insert an arbitrary-sized test particle into any 
two hard-sphere systems will be equal if the two distri- 
butions have the same first four moments Jl^] . In this 
sense the two systems may be termed "equivalent" . Since 
a binary mixture can always be chosen so as to match 
any four moments we look at the "equivalent" binary 
substitutionally-disordered crystal, for which simulation 
data is available |18|. By looking at test particles with 
sizes equal to the two species in the binary mixture, for 
which the chemical potentials are known, the remain- 
ing two unknown coefficients (Ai and A2) in the general 



expression for p ex (R) are determined. The resulting pre- 
dictions for the polydisperse crystal has been compared 
with simulation data previously Q . Agreement is good. 
For the polydisperse fluid accurate expression for 
p ex {R) are available. We use the approximate BMCSL 
]l9[ equation of state which for a Schultz distribution has 
the closed form 
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The excess free energy per particle 



is found by integrating Eq. |3[ Differentiation then yields 
an expression for the particle potential p ex (R) which is 
of the form of Eq. [| 

The total polydisperse free energy / (with / = F/V) 
consists of ideal and excess terms, / = f ld + f ex , which 
depend in a very different manner on the distribution 
p{R). The excess free energy, f ex = J dRp(R)p ex (R), is 
a function only of the four moments variables mo, . . . TO3. 
The ideal term f3f td = J dRp(R) \n(p(R)) , by contrast, 
depends upon the detailed shape of the function p(R) 
so formally, at least, the total free energy / resides in 
an infinite dimensional space. Sollich, Cates and Warren 
pO[ have shown that the full polydisperse phase diagram 
can be approximated by replacing the ideal free energy 
by a projected term / !d ({mJ) which includes only those 
contributions that depend on a finite set of moment vari- 
ables. The remaining contributions to the ideal free en- 
ergy, from those degrees of freedom of p{R) which can 
be varied without affecting the selected moments, are 
chosen to minimise the free energy. The power of this 
approach is that by including more moment variables 
the calculated phase diagram approaches, with increas- 
ing precision, the actual phase diagram. The position 
of equilibrium is fixed by the equality of the 'moment' 
chemical potentials, pi = df/dnii and the pressure P 
among all phases with / the projected free energy. For 
polydisperse hard spheres the excess moment chemical 
potentials are simply combinations of the (known) coef- 
ficients {A J in p ex {R) (Eq. |) since p{R) = Sf/6p{R) = 
J2i(df/dmi)wi(R) = y] f p jWj(R). The first two ideal 
moment potentials are [pO[ , ignoring unimportant fac- 
tors, p}f = lnp — ah\R and p l f = —aR. 

In order to understand the qualitative features of the 
polydisperse transition, we consider first the simplest de- 
scription in which only the lowest moment (too) is re- 
tained in the projected free energy. In this limit, there is 
no size fractionation so the normalised diameter distribu- 
tion, p(R), is fixed and equal in all phases. The location 
of the fluid-solid transition is determined by equating P 
and po, the chemical potential of the mean-sized particle, 
in each of the crystal and fluid phases. In this way we ob- 
tain the phase diagram of Fig. ftl At low densities we find, 
in qualitative agreement with previous work BH^,n2H14[ , 
that the density discontinuity at freezing Ap reduces with 



increasing polydispersity and eventually vanishes at the 
point at = 0.0833 and pt — 1.111. However, at high 
polydispersity, the calculated diagram contains a novel 
feature. For 0.07 < a < 0.083 we find a further transi- 
tion from the crystal back to a disordered phase M] . The 
location of this polydispersity- induced-melting transition 
varies sharply with polydispersity. The range of densities 
for which a crystal is found shrinks with increasing poly- 
dispersity until at at the crystal of density p t disappears 
completely from the phase diagram. At the point {pt^af) 
the line of fluid-to-crystal transitions intersects an upper 
line of crystal-to-amorphous transitions. At all points in 
the (p,a) plane the freezing transition remains first-order 
so the singularity at (pt,Q't) is equivalent to the point of 
equal concentration |22| seen in molecular mixtures and 
is not a critical point - so providing an answer to the 
second of our questions. 

We now turn to the vanishing density discontinuity 
in the vicinity of the point of equal concentration. The 
Gibbs free energy difference Ag — g s ~gi (with g = G/N) 
between the solid and liquid phases as a function of pres- 
sure for three fixed values of a is shown in Fig. g. The 
re-entrant nature of the freezing transition is very evident 
with a stable crystal appearing only in an intermediate 
range of pressures bounded by the two transitions where 
Ag = 0. The density change Ap at the liquid-solid tran- 
sition is given by the slope of the free energy curve at the 
point Ag = since dAg/dp = (l/p s ) — (l/l°0- Increas- 
ing the polydispersity raises the free energy of the solid 
relative to the fluid, displacing the Ag curve vertically 
and as is evident from Fig. [3 reduces the density jump at 
the transition. At the terminal polydispersity the solid 
just touches the fluid curve so the tangent is horizon- 
tal and Ap — 0. In a system of hard spheres (where 
the internal energy is constant) the condition Ap = 
necessarily requires the entropy change at this point to 
also vanish. Clearly while the underlying microscopic 
transition remains first-order the first derivatives of the 
thermodynamic potential are continuous at at- A con- 
ventional classification of this transition, following the 
ideas of Ehrenfest [g3|, is clearly inappropriate. 

Retaining more moments in the projected free energy 
allows the possibility of different-sized particles to be par- 
titioned between phases. To establish the effect of frac- 
tionation we have recalculated the phase equilibria with 
two moment variables. The phase diagram, now given 
by equating P and the moment potentials p$ and p\ in 
all phases, is unchanged in topology from Fig. [0. The 
point of equal concentration is retained although shifted 
slightly to (p tl a t ) = (1.115,0.0831). Hence our predic- 
tion of a re-entrant freezing transition seems to be robust. 
The extent of fractionation is generally small, although 
increasing as a — > a t , with the larger particles preferen- 
tially found in the crystal phase. Details of our calcula- 
tions are given elsewhere j24fl. 

The appearance of an equilibrium amorphous phase 



may be understood simply from maximum packing argu- 
ments. For uniform-sized spheres the maximum density 
of a randomly packed Bernal glass (p rC p — 1-22) is signif- 
icantly smaller than the geometric limit of a close-packed 
hexagonal of fee crystal (p cp = \/2). The greater pack- 
ing efficiency of the crystal ensures, that at high densi- 
ties, particles have more freedom and so a higher entropy 
than those in the fluid phase pa]. The stable high den- 
sity phase of uniform hard spheres is therefore crystalline. 
Polydispersity affects crystalline and disordered phases in 
different ways. In an amorphous phase, small particles 
pack in the cavities between large particles and p rcp in- 
creases with a while the constrained environment of a 
fixed repeating unit cell causes the maximum density of 
a crystal p cp to decrease with a. Computer simulations 
||ll| , P6| indicate that the limiting densities of amorphous 
and crystalline structures become equal at a « 0.05. 
For higher polydispersities disordered structures fill space 
more efficiently than ordered ones. Consequently the ap- 
pearance of an equilibrium amorphous phase and the en- 
suing re-entrant freezing transition should be a universal 
feature of all polydisperse systems - so answering the last 
of our questions. 

In conclusion we have presented a simple mean-field 
model of polydisperse hard spheres which suggests that 
the equilibrium state at high polydispersities and densi- 
ties is amorphous. An equilibrium crystal is found only 
at intermediate densities. The growing stability of the 
fluid phase with polydispersity causes a singularity in the 
density-polydispersity phase diagram which we identify 
as a point of equal concentration. Finally, although we 
use mean-field theory, our results should be robust with 
respect to fluctuation effects since the transition we find 
is not critical and the thermodynamic functions are not 
singular at this point. 
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FIG. 1. Phase diagram of a polydisperse system of hard 
spheres showing the re-entrant freezing transition. The den- 
sity discontinuity Ap = p s — pi vanishes at the point of equal 
concentration, marked by the filled circle. The inset figure 
shows the phase boundaries near this point in greater detail. 
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FIG. 2. The Gibbs free energy difference (per particle) Ag 
between crystal and fluid phases as a function of the pressure, 
for different polydispersities. The circles are the first-order 
fluid/solid transitions. The filled circle marks the point of 
equal concentration. 



